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Abstract 

Uptaking of contaminants by solid phases is relevant to many issues in environmental science as this process 
can remove them from solutions and retard their transport into the hydrosphere. Here we report on two structural 
studies performed on As-doped gypsum (CaSCU 2H2O) and calcite (CaCOs), using neutron (D20-ILL) and x- 
ray (ID11-ESRF) diffraction data and EXAFS (BM8-ESRF). The aim of this study is to determine whether As 
gets into the bulk of gypsum and calcite structures or is simply adsorbed on the surface. Different mechanisms of 
substitution are used as hypotheses. The combined Rietveld analysis of neutron and x-ray diffraction data shows an 
expansion of the unit cell volume proportional to the As concentration within the samples. DFT-based simulations 
confirm the increase of the unit cell volume proportional to the amount of carbonate or sulphate groups substituted. 
Interpolation of the experimental Rietveld data allows us to distinguish As substituted within the structure from 
that adsorbed on the surface of both minerals. Results obtained by EXAFS analysis from calcite samples show 
good agreement with the hypothesis of replacement of As into the C crystallographic site. 
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1. Introduction 

Arsenic is recognized as a dangerous pollutant of the 
environment [1] . Arsenic present in groundwaters may 
be trapped in the solid phase of minerals like calcite or 
gypsum, either by adsorption or by co-precipitation. 
When a contaminant is incorporated in the bulk rather 
than simply adsorbed at the surface, it is less available 
and it can be considered "immobilized" in the environ- 
ment at least until the host phase dissolution. The aim 
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of this study is to elucidate whether the incorporation 
of As(III) and As(V) into the bulk of calcite and gyp- 
sum, respectively, occurs or not, and to what extent. 

Surface-sensitive X-ray Standing Wave (XS W) stud- 
ies by Cheng et al. [2] show that the As atom re- 
places the C atom in the carbonate molecules of calcite. 
The geometry of the carbonate group is not preserved, 
showing a displacement of the As atom of 0.76 A in 
the [001] direction. Density Functional Theory (DFT) 
based simulations have proved that this replacement 
drives to a similar displacement of the As atom when 
C atoms are replaced by As atoms in the bulk of cal- 
cite (see below). This fact leads us to keep the same 
hypothesis of arsenite/carbonate replacement in our 
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study of As incorporation into the bulk. 

The stability diagram of aqueous solutions of 
H„As02~ 3 shows that the arsenate (AsO 3- ) is the 
most stable specie under oxidizing conditions. This 
anion is a tetrahedron with the As atom at the cen- 
tre, surrounded by four O atoms. The fact that both, 
arsenate and sulphate groups have the same geometry 
supports the hypothesis of a possible replacement of 
sulphate by arsenate groups when gypsum is precipi- 
tated in presence of As(V). The charge is compensated 
by bonding to an extra proton. 



2. Materials and methods 

In order to test the possible mechanisms for As im- 
mobilization by calcite and gypsum, samples of both 
minerals were synthesized in the presence of As(III) 
and As(V), respectively. Calcite precipitation was con- 
ducted at pH = 7.5 by addition of CaCL; and Na2C03 
solutions. Gypsum was precipitated from Na2S04 and 
CaCl2 solutions at three different pH values: 4, 7.5 and 
9. Arsenic concentrations incorporated in the solids 
range between 30 mM/kg and 1200 mM/kg for calcite 
and 100 mM/kg to 1000 mM/kg for gypsum. 

Powder samples were analysed by neutron diffrac- 
tion at the high flux powder diffractometer D20 
(ILL). Experiments were carried out using a Cu(200) 
monochromator which gives a wavelength of A = 1.30 
A and at ambient conditions of pressure and tem- 
perature [4]. Diffraction patterns were taken for the 
samples in their container and for the empty cell, in 
the range of 10° to 130°. Also powder diffraction ex- 
periments were performed with x-ray at ID11 (ESRF), 
reproducing the same experimental conditions but us- 
ing a wavelength of A = 0.5 A. Both diffraction data 
sets were analysed using Fullprof [5]. 

Geometrical optimisations of the unit cell and the 
supercells of pure and As-doped calcite and gypsum 
were done with the Vienna Ab-initio Simulation Pack- 
age (VASP) [6]. The PBE functional and PAW pseu- 
dopotentials were used. The goal was to reproduce the 
expansion of the unit cell produced by the incorpora- 
tion of As atoms into the structure of both minerals. 
Unit cells of pure calcite and gypsum obtained from 
Rietveld refinements were used as starting point for all 
the models. Geometrical optimizations of single unit 
cells and of supercells of gypsum were done replac- 
ing the sulphate molecules SO4 - by arsenates AsO 3- . 
Similar simulations were performed with supercells of 
calcite replacing C by As atoms [4] . 

EXAFS data were collected on a diffraction and 
absorption beamline (GILDA-BM8) at the ESRF of 
Grenoble and extracted using standard procedures [7]. 
The theoretical photoelectron paths were generated 
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Fig. 1. Experimental unit cell volume of As-doped calcite and 
gypsum samples in function of As concentration. 

using the FEFF8 [8] code and the fit performed using 
the MINUIT library from CERN [9]. A monochroma- 
tor of Si(311) was used to set the incident energy at 
the K-edge of As (11867 eV). 



3. Results and discussion 

Our diffraction data show an expansion of the unit 
cell due to As incorporation into calcite crystallites 
(Fig. 1). By modelling, As concentration in the bulk 
of the samples can be extrapolated by comparison of 
the relative volume changes between the experimental 
and simulated data. Calcite unit cell and two super- 
cells (2x2x1 and 3x2x1) were geometrically optimised 
replacing one and two units of AsO^ - (150 mM/kg and 
290 mM/kg) by COg~ units. The simulations showed 
a volume expansion linearly dependent on the replace- 
ment of As, as expected by Vegard's law (Fig. 2). This 
augmentation is due to the lattice expansion along the 
a axis as As concentration in solids increases. The ob- 
served As ion displacement of 0.57 A over the O base 
along the (0001) is compatible with Cheng's results [2]. 
The lower value of this displacement (0.57 A vs. 0.76 
A) can be due to the fact that atoms near the surface 
are less attached to the solid and can move more freely. 

The experimental value of the unit cell volume was 
interpolated using the linear fit of the simulated vol- 
ume expansion, giving values of 9, 10 and 16 mM/kg of 
As in calcite for the three synthesised samples. Simula- 
tions of 3x2x1 supercells with one CO3 - unit replaced 
by one AsOg - unit were done to check whether the re- 
placement is more likely in sites at the same or at dif- 
ferent crystallographic planes. The higher enthalpy of 
formation (AH — 211.36 meV) for the calcite struc- 
ture with two As atoms lying on the same plane shows 
that replacement is more likely to happen in different 
planes, leading to a more stable structure. 
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Fig. 2. Simulated relative expansion of unit cell volume of gyp- 
sum and calcite. The insets show a detail of the interpolation 
of the experimental data using the simulated expansion. 



We found from EXAFS data analysis a nearest 
neighbour distance of w 1.77(3) A corresponding 
to the As-O bond distance. This value lies in be- 
tween the reported one [10] for the arsenite molecule 
(dAs-o = 1-84 A) and that obtained from the simula- 



tions (d%™ — 1.75 A). The coordination number was 
kept fixed to its theoretical value (TV = 3) in order to 
reduce the correlation between free parameters in the 
fitting procedure. This result supports the hypothe- 
sis of the incorporation of the As atoms into the C 
crystallographic sites. 

Figure 1 shows the gypsum unit cell volume obtained 
from combined refinements of neutron and x-ray data. 
The expansion of the unit cell is proportional to the 
As concentration in solids and strongly dependent on 
the pH value: the biggest expansion is found in sam- 
ples synthesized at pH 9. This result is in good agree- 
ment with the hypothesis of replacement of sulphate 
(SO^") by arsenate groups (AsO^ - ). This replacement 
is more likely at higher pH values, according to avail- 
able thermodynamical data regarding speciation of As 
[3] . The expansion of the unit cell parameters is due to 
the different lengths for As-O (« 1.70 A) [12] and S-O 
(« 1.50 A) [11] bonds. 

Simulations show an increasing of the unit cell vol- 
ume proportional to the number of atoms of S replaced 
by As (Fig. 2). Four single cells (with 940, 1809, 3357 
and 4696 mM/kg of As) and four supercells were simu- 



lated: two 2x1x2 supercells with one and two As atoms, 
a 2x1x3 and a 3x1x3 with one As atom each, giving 
As concentrations of 358, 705, 240 and 160 mM/kg, re- 
spectively. The simulations allow us to extrapolate the 
As concentration in the bulk of the samples by com- 
paring the relative volume variations between the ex- 
perimental and simulated data (Table 1). 
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Table 1 

Extrapolated values of the As concentration in samples of gyp- 
sum. The result for pH 4 is due to the near zero relative ex- 
pansion of the unit cell. 

Our results support the hypothesis of As immobili- 
sation by incorporation into the bulk of these minerals. 
This improves the knowledge on the long term stabil- 
ity of contaminated sludges and it has important con- 
sequences for site remediation actions. This work is an 
example of a direct link between fundamental research 
and environmental issues. The understanding of the As 
compounds behaviour in sedimentary environments is 
essential to estimate and predict possible consequences 
of forecast or accidental events. 
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